from Excitation_Carbon import *

cs = CS_Excitation(species_name = "C+3")


id = 0

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =8.0,
    A   = 4.744,
    B   = 4.440,
    C   = 5.971e-1,
    D   = 0.0,
    E   = 4.519,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+3_2s2S-2p2P")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =37.6,
    A   = 4.846e-1,
    B   = -2.743e-1,
    C   = 4.768e-1,
    D   = -3.971e-1,
    E   = 0.0,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+3_2s2S-3s2S",
    X1  = 1.072,
    P   = 7.390e-1,
    Q   = 2.460e-1)


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 40.0,
    A   = -5.731e-1,
    B   = 9.548e-1,
    C   = -8.931e-2,
    D   = 0.0,
    E   = 5.638e-1,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+3_2s2S-3p2P",
    X1  = 1.019,
    P   = -3.412,
    Q   = 9.141)


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =40.0,
    A   = 1.252,
    B   = -1.243,
    C   = 4.93e-1,
    D   = 1.326e-1,
    E   = 0.0,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+3_2s2S-3d2D")


id = id + 1
cs.add_cross_section(
    id = id, 
    Vif =50.6,
    A   = 2.775e-1,
    B   = -6.663e-1,
    C   = 4.959e-1,
    D   = 0.0,
    E   = 3.530e-6,
    F   = 0.0,
    omega = 2.0,
    equation_number = 6,
    reaction_name = "Excitation_C+3_2s2S-4p2P")



cs.plot_cs()
cs.save_cs()
cs.cal_total_cs()